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ABSTRACT 

We study atomic line diagnostics of the inner regions of protoplanetary disks 
with our model of X-ray irradiated disk atmospheres which was previously used to 
predict observable levels of the Nell and NeHI fine-structure transitions at 12.81 
and 15.55 /im. We extend the X-ray ionization theory to sulfur and calculate the 
fraction of sulfur in S, S + , S 2+ and sulfur molecules. For the D'Alessio generic T 
Tauri star disk, we find that the SI fine-structure line at 25.55 /im is below the 
detection level of the Spitzer Infrared Spectrometer (IRS), in large part due to 
X-ray ionization of atomic S at the top of the atmosphere and to its incorporation 
into molecules close to the mid-plane. We predict that observable fluxes of the 
SII 6718/6732 A forbidden transitions are produced in the upper atmosphere at 
somewhat shallower depths and smaller radii than the neon fine-structure lines. 
This and other forbidden line transitions, such as the O I 6300 / 6363 A and the 
CI 9826/9852 A lines, serve as complementary diagnostics of X-ray irradiated 
disk atmospheres. We have also analyzed the potential role of the low-excitation 
fine-structure lines of C I, C II, and O I, which should be observable by SOFIA 
and Herschel. 

Subject headings: accretion, accretion disks - infrared: stars - planetary systems: 
protoplanetary disks - stars: formation - stars: pre-man- sequence - X-rays: stars 
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Introduction 



There is growing evidence for the existence of warm atmospheres in the inner disks of 
low and intermediate mass young stellar objects (YSOs). They first became apparent from 
observations of the CO overto ne bandhead e missio n near 2.3 /im in Herbig Ae stars and 



T Tauri stars (e.g., Carr Il989l ; Najita et al. |2000| ). More recently, the emission from the 



fundamental CO ro - vibrational band near 4.6 /xm has been observed in many T Tauri stars 
(|Naiita et al.l 120031 ) . as has the UV fluorescence of molecular hydrogen (e.g., Herczeg et al. 
20021 ). These observations indicate that temperatures in the range 1000-2000 K occur over 



a considerable vertical column density of hydrogen ~ 10 cm in the inner disk region. 
This and other spe ctroscopic evidence for warm gaseous disk atmospheres were reviewed by 
Najita et al.l (120071 ). where extensive references are given to earlier work. 



Calvet et al.l (Il99ll ) pioneered the now widely accepted view that disk atmos pheres 



are heated by stellar radiation. In previous publications (Glassgold, Najita & Igea 12004 
henceforth GNI04) we argued that stellar X-rays play an important role in heating the 
upper gas layers of circumstellar disks on the following basis. First, significant X-ray emission 
appears to be a universal aspect of low-mass YSOs; second, X-rays couple strongly to the gas; 
and third, hard X-rays with energies > 1 keV can penetrate deeply into disk atmospheres. 

One special aspect of X-ray irradiation is that temperatures as large as 4000-5000 K are 
achieved high in the atmosphere. We recently pointed out that this hot layer gives rise to 
signifi cant flux levels of the fine-structure lines of neon at 12.81 /j,m and 15.55 /zm (Glassgold 
et al. 120071 . henceforth GNI07). 



The Nell 12.81 /im line has now been detecte d in many YSOs by t 



Spectrometer (IRS) in nearby star-forming regions (IPascucci et al. 



2007 



le Spitzer Infrared 



Lahuis et al. 2007 



Espaillat et all 120071 ) and by MICHELLE at Gemini North flHerczeg et al.l 120071 1 . The mea- 
sured fluxes are in good agreement with the predictions in GNI07. These observations 
support the hypothesis that disk atmospheres are heated and ionized by photons energetic 
enough to ionize Ne (which has first and second ionization potentials of 21.56 and 41.0 eV), 
most likely X-rays. A key motivation in the work of GNI07 was the search for clear diagnos- 
tics of the effects of X-ray irradiation. The Ne lines seemed particularly appropriate because 
hard X-rays generate high ionization states of neon via the Auger effect and because neon 
chemistry is especially simple. Gorti & Hollenbach (2006, private communication) have 
suggested that EUV radiation may contribute to the neon fluxes, although the emission 
characteristics of YSOs in this wavelength band are poorly known compared to X-rays. 



In this paper we extend the considerations of GNI07 and seek other fine-structure and 
forbidden transitions that have a potential to complement the neon lines. We have been 
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guided in part by the intrinsic strength of a line, which we can take as the thermalized 
emissivity in a transition from an upper to a lower level u — > I, 

j(u, I) = P{u) xn n A{u, I) E{u, I). (1) 

Here inn is the density of the emitting species, P{u) is the population of the upper level, and 
A(u, I) and E(u, I) are the Einstein A-value and energy of the transition. The last two factors 
in equation ([T]) emphasize high frequency transitions, but the other factors are crucial since 
they require a theory for the ionization of the carrier species (the x factor) and the excitation 
of the level (the factor P(u)), since they are sensitive to the physical conditions. We have 
chosen to focus here on the lines of oxygen, sulfur and, to a lesser extent, carbon. Some 
spectral lines of particular interest are the 25.25 /zm fine-structure transition of SI, which 
has been searched for, mostly without success, by Spitzer, the fine-structure lines of O I, 
which should be detectable with SOFIA and the PACS imaging spectrometer on Herschel; 
and various optical-infrared forbidden transitions accessible from the ground. 

The atomic carriers of these lines appear in the many disk chemical models developed 
over the last 10-15 years (reviewed by Bergin et al. 120071 ). They are closely related to the 
photon-dominated models (PDR) developed for the interstellar medium. In this case stellar 
and interstellar radiation generate warm atomic regions in the upper and lower layers of the 
disks until the UV radiation is absorbed and molecule formation becomes efficient. Jonkheid 



and Kamp and their collaborators (IJonkheid et al.ll2004l . 120061 120071 : iKamp et al.ll2003l ) have 
gone further and calculated the emissivity and lineshape of the fine-structure lines of C I, 
CII, and OI. 

The 63 fim fine-struc ture l ine of O I was detected in YSOs w ith th e KAO (e.g., Co hen et 



al. Il988l ; Ceccarelli et al. 
Creech- Eakman et al 



2002 



1993) and with IS O (e.g ., Nisini et al. Il999l : Spinoglio et al. 12000 



Liseau et al. 120061 ). but the interpretation of the results is 
hampered by the low spatial and spectral resolution of the observations so that it is unclear 
what fraction of the emission arises from the disk. Disks have also been considered as 



a possible source of forbidd en line emission by lHartigan et al.l ( 119951 ). iKwanl (119971 ) and 
Storzer fc Hollenbachl (120001) . The forbidden lines have been well studied in YSOs (see, e.g., 
the review by Ray et al. 120071 ). and some of the observed "low- velocity-compo nents" of th e 
emission may originate from the disks in these systems (e.g., Kwan Sz Tademaru Il988l . Il995l ) . 
We will return to the discussion of the observations of the fine-structure and forbidden lines 
in Section 5. 

The calculations reported in this paper are exploratory in nature, as were those in 
GNI07. They examine the ability of X-rays, a well characterized property of TTS, to produce 
atomic line emission from T Tauri disks. The model employs the continuous disk density 
distribution developed by D'Alessio et al. (1999) for a generic T Tauri star, and they feature 
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stellar X-ray heating and ionization in the context of a simplified thermal-chemical model, 
which will be discussed in more detail in Section 2. We do not attempt to fit the observations 
for any specific system, but instead try to derive conclusions that might be applicable to any 
protoplanetary disk. 

The plan of the rest of this paper is as follows. In the next section, we review and 
update the thermal-chemical model of GNI04. In Section 3, we summarize the essentials of 
the basic ionization theory and extend it to include the X-ray ionization of sulfur. In Section 
4, we outline the excitation and flux calculations for the lines of interest, and then discuss 
the relevant observations in Section 5. The last section summarizes the main results of this 
paper. 



Review Of The Model 



We use the thermal-chemical model of GNI04 with minor corrections and updates. The 
important processes in determining the thermal balance are X-ray heating, gas-grain heating 
and cooling, mechanical heating (e.g., from disk accretion or the wind-disk interaction), 
and line cooling (Lya, recombination lines, 01 fine-structure and forbidden lines, and CO 
rotational and ro- vibrational lines). The chemical network contains about 125 reactions and 
25 species. Errors in the specification of these processes have been corrected, and the entire 
reaction base has been re-evaluated and updated. None of the changes affects the conclusions 
of GNI04 and GNI07 in any major way. GNI04 provide a detailed description of the physical 
processes, and the chemical reaction base is available from the authors on request or can be 



downloaded at http: / / astro.berkeley.edu/~rowin/MGN_rates.pdf| 



Following GNI04 and GNI07, we calculate the chemical and thermal structure of the 
gas in a protoplanetary d i sk for the smooth density distribution of the generic T-Tauri disk 
model oflD'Alessio et al.l (119991 ) using an updated version of the code developed by GNI04. 
D'Alessio et al.l (Il999l ) did not distinguish between the thermal properties of the dust and 
the gas, and their density structure is in hydrostatic equilibrium. However, at high altitudes 
and low vertical column densities, the temperature of the gas is much higher than the dust 
(GNI04), and this implies that our model is no longer in hydrostatic equilibrium. Were we to 
treat the gas and dust as two separate fluids and require hydrostatic equilibrium, the results 
would differ somewhat from the exploratory calculations reported here and in GNI07. For 
example, we would expect the gaseous disk to be puffed up close to the star and to shield 
the outer part of the disk from stellar radiation. On the other hand, the flaring of the outer 
disk would also be increased , and thus it would be illum inated more directly by the central 
star (jPullemond et al.l 120021 ; iGorti fe Hollenbachl 120041 ). A quantitative evaluation of such 
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effects on the diagnostic line ffuxes must await the development of a more complete model 
that treats both the gas and the dust self-consistently. 



The D'Alessio model does not include X-rays. In our previous studies based on this 
density distribution, we adopted a "reference" X-ray luminosity of L x = 2 x 10 30 erg s _1 
and a thermal spectrum with kTx = 1 keV. It is well known that both the X-ray luminosity 
and spectrum of young stellar objects can change with time, and that they vary widely from 
source to source. The i onization rate is no t very sensitive to the temperature of the spectrum, 
as shown in Fig. 5 of llgea fc Glassgoldl (119991). In the Chandra O rion Ultradeep Project 
(COUP) study of the Orion Nebula Cluster (IGetman et al.l 120051 ). the X-ray luminosity 
function for the entire cluster spans 5 decades from \ogL x = 27 — 32. Two of these decades 
of variation may be associat ed roughly with the d ependence of Lx on stellar mass and one on 
a dependence on stellar age (IFeigelson et al.ll2005l ). leaving two or more unexplained decades. 
Thus the X-ray luminosity of a pre-main-sequence object with specific mass and age may 
vary over two orders of magnitude. This is consistent with the findings of I Wolk et al.l (120051 ) 
for a subset of 28 solar-mass stars in the Orion Nebula Cluster with masses in the range 
M* = 0.9 — 1.2M Q , which have a median X-ray luminosity of \ogLx = 30.25 erg s _1 . The 
spread of two orders of magnitudes can only be partly explained by variability, including 
flares. Similar results hay e been found in the X MM- Newton study of the Taurus- Aurigae 
cluster ( IGiidel et al.l 120071 : iTelleschi et al.l 120071 ) , where the X-ray luminosity ranges over 
three decades (logLx = 28 — 31). 

Our previous choice Lx = 2 x 10 30 erg s _1 was consistent with the observations of 
young T Tauri stars with masses ~ j-M^ . an d with the X-ray lumin osity of well studied 
sources such as T W Hya (IKastner et al.l 120021 ; IStelzer fc Schmittl 12004 ) . The median X-ray 
luminosity measured by XEST for classical T Tauri stars (with disks) in the Taurus- Aurigae 
cluster is Lx = 5 x 10 29 erg s -1 , four times smaller than our reference value. In recognition 
of this difference as well as in the observed 1-2 dex spread in measured luminosity, we will 
consider a range of X-ray luminosities between Lx = 2 x 10 29 — 2 x 10 31 erg s" 1 . X-ray 
luminosities in this range are likely to yield fluxes for the lines of interest that can be 
detected with current facilities. 

In order to improve the numerical accuracy of the flux calculations in Section 4, the 
calculations reported here cover an extended range of radii from 0.25-100 AU, in contrast 
to GNI07 who considered the range 1-40 AU. In making this extension, we do not consider 
any structural and physical modifications that might occur at very small and very large 
radii. For example, spectral energy distributions measured with the Spitze r IRS indicate the 



occur rence of inner holes and rims at distances of the order of several AU (iDullemond et al. 



20071 ). Since we use the D'Alessio model calculations, which varies continuously from 0.028 



- 6- 



to > 500 AU, the present calculations do not take into account structures such as holes, gaps, 
and rims. 

The elemental abundances used here and in the flux calculations to follow are: x^ e = 0.1, 



x c = 2.8 x 10~ 4 , x = 6.0 x 10" 



7.0 x 10" 6 , 



1.0 x 10~ 4 , and x Na = 1.0 x 10~ 6 . 



The abundance of sulfur that we have adopted is in the range suggested by the recen t 
modeling of sulfur molecules observed in the Horsehead Ne bula (Goicoechea et al.ll2006l ). 
It is much closer to the solar photospheric/meteoritic value (lAsplund et al.l 120051) th an the 
heavily depleted values used in modeling dark clouds (e.g., Millar & Herbst Il990l ). UV 
absorption line studies of the diffuse interstellar medium provide no eviden ce for depletion of 
sulfur, as they do for more refractory elements (jSavage Sembachlll996l ). Fig. [1] shows the 
thermal structure of the reference disk model. The stellar and disk parameters that define 
the D'Alessio 1999 model are given in the caption. For high altitudes (Ah < 10 21 cm -2 ) 
and moderate radii (R < 25 AU), temperatures as high as 4000-5000 K are reached. Going 
deeper into the atmosphere, T undergoes a sharp drop toward mid-plane temperatures in 
the range from 30 to 200 K, depending on radius. Near the mid-plane, where the density 
is high (n = 10 8 — 10 10 cm -3 ), the gas and the dust are thermally coupled. For these radii 
(< 25 AU), the temperature starts high at the top of the atmosphere and then increases 
further with increasing vertical column density or decreasing altitude (see also Figure 2 of 
GNI04). This is a consequence of the balance between X-ray heating and Ly-a cooling 
that controls the temperature at the very top of the atmosphere. The heating and the 
cooling are both proportional to the first power of the volume density. The small increase 
in temperature arises from its logarithmic dependence on optical depth, which of course 
increases with increasing vertical column density. Beyond 25 AU, the O I fine-structure lines 
are the dominant coolants at the disk surface. 



3. Ionization theory 

The GNI07 thermal-chemical program includes the elements H, He, C and O, plus a 
generic heavy atom represented by Na. X-ray ionization of H, He, H2 and C are included 
explicitly, and that of O implicitly. Generally speaking, X-ray ionization of a cosmic gas 
occurs by the absorption by the K and L shell electrons of heavy atoms and ions. The 
resultant photo and Auger electrons then generate many more secondary electrons by the 
collisional ionization of H and He. The ionization state of a heavy atom A or ion in a mainly 
atomic gas is determined by several processes: direct X-ray ionization (rate Cdir); secondary 
electron ionization (rate Csec); electronic recombination; and charge transfer to H and He. 
The latter process may be fast (rate coeffcient ~ 10 -9 cm 3 s _1 ) or slow. In the case of neon 
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(GNI07), charge transfer is relatively slow and fosters high abundances of neon ions. The 
ionization rates for atom A, dir(A) and ( sec (A), are always per atom, whereas the total 
ionization rate ( is per H nucleus. 

On the basis of its elemental X-ray absorption cross section, we estimate that the oxygen 
ionization rates are Cdir(O) = 17( and ( sec (0) = 2A(, where ( is the X-ray ionization rate of 
the cosmic mix of disk gas. However, near-resonant charge exchange of H + and O is fast and 
dominates the ionization of oxygen. Thus direct X-ray ionization of oxygen can be ignored, 
as in GNI04. In the case of carbon, we estimated Cdir(C) = 6( and Csec(C) = 4£, but now 
charge exchange of H + and He + with C is very slow, and therefore X-ray ionization of carbon 
is important and is included explicitly in the thermal-chemical program. 

Figure [2] shows the electron fraction calculated for our model; it exceeds 0.01 at the top 
of the atmosphere down to a depth of ~ 10 19 cm~ 2 . The electron fraction decreases more 
smoothly with vertical column than the temperature. Near the mid-plane, its value is larger 
than that given by GNI07 because io nization due to the radioactive decay of 26 Al have bee n 



included atarateC 26 = 4xl0~ 19 s" 1 flStepinskilll992l : lGlassgoldlll995l : lFinocchi fc Gaillll9971 ). 



This assumes that the entire disk is optically thick with respect to the 1.809 MeV decay 
7-rays, and that the 26 A1/ 27 A1 ratio has the canonical value found in meteoritic calcium- 
aluminum inclusions (5 x 10~ 5 ). It also ignores the effects of the lMyr mean life of 26 Al 
and the possibility that the initial distribution of 26 Al is spatially inhomogeneous. In light 
of these simplifying assumptions, it is likely that the rate is even smaller than the one 
we have adopted, but the exact value is not crucial for the main subject of this work^. 
We use a separate program to obtain the abundances of neon and sulfur species. The 
ionization calculation of neon is essentially the same as GNI07. Figure [3] shows the fractional 
abundances of Ne + and Ne 2+ relative to the total abundance of neon. Abundances of the 
order ~ 0.1 or more are found at high altitudes (Nh < 10 21 cm -2 ). It is remarkable that 
X-ray ionization produces significant abundances of Ne + and Ne 2+ out to radii as large as 
30 AU. 

The closed shell character of neon means that the complexities of molecule formation 
and destruction do not enter, and only three species, Ne, Ne + and Ne 2+ have to be considered 
in estimating their line fluxes. This is not the case for the calculation of the S I and S II 
fine-structure and forbidden lines. In order to calculate the abundances of S and S + , we 
have to consider the transition into sulfur molecules, such as SO, SO2 and CS, which occurs 
at large vertical column densities. Figure H] gives a schematic overview of sulfur ionization 



1 Many modelers of disk chemistry have used 26 Al ioniz ation rates an order of magnit ude larger than we 
have as the result of a typographical error in the paper by Umebavashi fc Nakanol ( 1981 ). 
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including a simplified warm chemistry, following iLeen fc Grafa (119881 ) 



Sulfur ions in the disk are mainly produced by X-ray ionization. Sulfur has L and K 
edges near 0.2 and 2.5 keV and, since we consider X-ray spectra that are not particularly hard 
(Tx ^ 1 keV), photons with energies between the L and K edges are the ones most strongly 
absorbed by sulfur. Because charge exchange with H is fast for highly-ionized sulfur ions, 
X-ray ionization of sulfur accompan i ed by the A uger process leads prim a rily to S 2+ and S 3+ . 
Calculations by iButler fc Dalgarnol (119801 ) and IChristensen fc Watsonl (119811 ) suggest that 
S 3+ charge exchange with atomic hydrogen is also fast, while that of S 2+ is slow. Therefore, 
we only include ions up to S 2+ , as we did in the case of neon (GNI07). However, in addition 
to electronic recombination and charge exchange of sulfur ions with atomic hydrogen, we also 
include the S + H + charge-exchange reaction, which can ionize S at a si gnificant rate below 
10000 K according to the theoretical calculations of IZhao et al. d2Q05[) Furthermore, the 
reaction S 2+ + H2, which has been measured to be fast (IChen et al. I l2003h . can significantly 
reduce the abundance of S 2+ . 

As for carbon and oxygen, we can estimate the direct X-ray ionization rate of sulfur 
from its X-ray absorption cross section. Its energy dependence is very similar to that of the 
mean X-ray absorption cross section averaged over interstellar or solar abundances. Thus the 
direct ionization rate of sulfur, resulting from the (L to K) 0.2-2.4 keV band is basically given 
by the ratio of the two absorption cross sections at these energies. For example, at 1 keV, the 
ratio is ~ 500. This then is the approximate ratio of the primary or direct X-ray ionization of 
sulfur (per sulfur nucleus) to that of the mean interstellar atom (per hydrogen nucleus, since 
abundances are conventionally normed to hydrogen). However, most of the ionization of the 
gas as a whole comes from secondary-electron ionization of hydrogen (atomic and molecular) 
and helium. Since there are roughly 25 such electrons produced per primary ionization, the 
ratio of the direct sulfur ionization rate per sulfur nucleus to the total X-ray ionization per 
H nucleus is 500/25= 20: 



C(S) dir ~ 20c. 



(2) 



The secondary electrons generated by X-ray ionization of all of the cosmic elements 
collisionally excite and ionize S at a rate that is roughl y given by the ratio of the electronic 
ionization cross section of sulfur to that of hydrogen (IMaloney et al.lll996l ). This leads to 
the approximate value, 

C(S)sec ^ 50 (3) 

The direct ionization rate is mainly responsible for producing higher S ions, which we lump 
together into S 2+ because of fast charge transfer. The secondary ionization rate primarily 
produces an additional electron, e.g., it generates S + from S. Some of the sulfur ion rate 
coefficients are listed in Table [TJ Those for radiative recombination are well calculated. As 
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in the case of neon, however, the rate coefficients for electron transfer from atomic hydrogen 
are not well established. They represent an important uncertainty in the results of this 
paper. 

To address the question of whether S remains the dominant species near the mid-plane, 
we include the formation and destruction of the sulfur molecule s, SO, SO2, CS and OCS, by 



using the warm neutral sulfur chemistry of (ILeen fc Grafilll988l ). We also include X-ray ion 



ization and destruction of these molecules, assuming that X-ray cross sections for molecules 
are additive over the constituent atoms, and using the atomic rates discussed above. We also 
make the simplifying assumptions that a direct absorption of an X-ray by a sulfur molecule 
always leads to dissociation, and that secondary electronic ionization always produces a 
molecular ion. The sulfur molecular ions are rapidly destroyed by dissocia t ive re combina- 



tion. The branching ratios are summarized by iFlorescu-Mitchell fc Mitchell! (120061 ) 



Figure [5] shows the fractional abundances of the sulfur species S + , S, SO2 and CS as a 
function of vertical column density above the disk mid-plane at a radial distance of 20 AU 
from the star. We find that sulfur is completely locked up in molecules for column densities 
larger than 0.5 — 1.0 x 10 22 cm -2 . In this region, neutral formation reactions start to dominate 
over X-ray destruction as the X-rays get shielded and the increased density enhances the 
reaction rates. The main species at large column densities is SO2, but there may also be a 
layer of CS at intermediate column densities. 

Figure M shows the S and S + fractions relative to the total abundance of sulfur in the 
disk for radii between 0.25 and 40 AU. S + is the dominant species at the highest altitudes 
for R < 25 AU , whereas S assumes that role at intermediate altitudes. Beyond R > 25 AU, 
X-ray ionization of sulfur is not strong enough to maintain S + as the dominant species. At 
low altitudes, we find a clear cut-off in the atomic S abundance, with almost all of the sulfur 
in molecules. 



4. Excitation and Flux Calculations 

The flux calculations for the atomic ions considered in this paper are relatively straight- 
forward because one needs to deal only with a few levels of relatively low excitation, even for 
the forbidden lines. For neon, we continue to focus on the ground fine-structure transitions 
of Ne + and Ne 2+ . For O, S and C, we treat the 5 lowest levels in the ground electronic 
configuration: a fine-structure triplet, that produces mid and far IR lines, and two higher- 
energy singlet levels that give rise to forbidden optical-IR transitions. For C + , we consider 
only the ground level fine-structure doublet that generates the famous 158 /im line. Simpli- 
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fications often occur. For example, the fine-structure transitions of CI and CII have low 
critical densities and, because of the high disk densities, are near-thermalized. Similarly, 
the forbidden lines are usually optically thin and the effects of line trapping are small. In 
the following sections, we discuss specific issues that arise in the calculation of the emis- 
sivities of the lines of interest. We use the line frequencies and A- values in the NIST data 
base ( ]http:/ /physics.nist.gov/PhysRefData/A SD/ind ex.html[ ). Electronic excitation rates 
are now reasonably well established, but those for excitation by H, H 2 , He and H + are often 
unknown, except perhaps for O I, C I and C II. References for our choice of rate coefficients 
are given in Table [2j In those cases where trapping plays a role, e.g., for the 01 and SI 
fine-structure lines, our results apply to a face-on disk. 



4.1. Neon 



We calculate the optically-thin Nell and Nelll lines fine-structure lines at 12.81 and 
15.55 /im in the manner of GNI07 with the following change. While retaining the two-level 
population formula for the ground- level doublet of Nell, we solve exactly for the populations 
of the ground-level triplet of Ne IIJ0. This more complete excitation calculation leads to an 
increase in the population of the upper J = level and, to a lesser extent, the J = 1 level, 
and increases the line intensities. 

Figure [7] shows the distribution of the Nell 12.81 /im emissivity throughout the disk 
for the reference X-ray luminosity. Most of the emission is produced at high temperature 
(~ 4000 K) and high electron abundance, corresponding to radial distances R < 25 AU and 
vertical column densities Nh < 10 21 cm™ 2 . The largest emissivity does not always occur 
at the highest altitude above the mid-plane. This is due to an increase of the Ne + density 
and its specific emissivity with increasing depth. The upper J = 1/2 level is sub-thermally 
populated and its population is proportional to the electron density, which also increases 
going down into the atmosphere. 

The above calculations assume that the excitation is solely due to electronic collisions. 
However, we have been reminded that, as the electron fraction decreases with height, colli- 
sions with atomic hydrogen begin to play a role (D. J. Hollenbach 2006, private communi- 



2 GNI07 omitted the weak quadrupolar transition that connects the top (J = 0) and bottom (J = 2) 
levels. Their approximate equation (3-4) also contained a typo; the correct formula is, 

p (o) = — (A) 

uy ' l + 5/3Ci_ 2 exp925.3/r+l/3exp(-399/T)/C _i v ; 
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cation). Using the results of iBahcall fc Woll (119681 ). Irlollenbach fc McKed (Il989l ) estimated 
the rate coefficient for the de-excitation of the 12.81 /im transition in collisions with atomic 
hydrogen to be 1.3 x 10~ 9 cm 3 s -1 . Although the theory of Bahcall and Wolf may give the 
right order of magnitude (as it does for the excitation of the O I fine-structure transitions), 
it is incorrect in principle since it implies that the excitation cross section is proportional to 
the scattering cross section. According to Bahcall and Wolf, the rate coefficient in a collision 
with a neutral atom varies with temperature as T 1//6 (c.f. the van der Waals potential) and 
as a constant in a collision with an ion (c.f. the Langevin potential). However, collisional 
excitation of a fine-structure transition requires an exchange of a spin between the incident H 
atom and the target electrons and thus depends on differences between potentials, differences 
which decrease more rapidly than 1/r 6 and 1/r 4 for neutral and ionic collisions, respectively. 



The full quantum-m echanical calc ulations for 01 and C I (ILaunav fc Rouefll Il977al ; 
Abrahamsson et al. 2007 ) and for CII ( Launay fc Rouefi 1977b ; Barinovs et al. 20051 ) bear 



out the fact that the Bahcall and Wolf theory predicts the wrong temperature dependence 
for the collisional exc itation rate coe f ficient s by atomic hydrogen. We have compared the 
recent calculations of iBarinovs et al.l (120051 ) for H + C II with the Bahcall and Wolf the- 
ory, and find that the latter over-estimates the rate coefficient by factors of 7.5 to 5 in the 
temperature range from 500-4000 K. If something like this reduction factor also applies to 
H + Nell, then the order of magnitude of its rate coefficient for collisional de-excitation 
of the Nell fine-structure doublet is ~ 2 x 10 cm 3 s . We have calculated the 12.81 /xm 
emissivity (including line trapping) with this value for the H de-excitation rate, and find that 
it increases the flux by about a factor of two. This increase arises mainly at large radial dis- 
tances, where the deviations from a thermal population are the greatest. We conclude that 
our calculations based on electronic excitation alone underestimate the Nell 12.81 /im flux, 
but possibly not by a large factor. Actual calculations of the atomic hydrogen excitation 
rate would be most welcome. 



4.2. Oxygen 

Figure [8] gives the 5 levels that we consider in calculating fine-structure and forbidden 
line emission from atomic oxygen and sulfur. According to Table [2], collisional excitation 
rates for the fine- structure levels of O I are available for collision partners e, p, H, H 2 and He. 
Atomic hydrogen tends to dominate and, even without line trapping, the critical densities 
are modest, typically n cr (H) < 10 6 cm -3 . At radii (> 25 AU), densities of this order are not 
reached until vertical column densities N-r ~ 10 20 — 10 21 cm -2 . The full emissivity calculation 
is particularly relevant at large radii. 
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Figure [9] shows the spatial distribution of the emissivity of the O I 63 fim fine-structure 
line. The results for the O I 146 fim line (not shown) are similar. The spatial profile is 
particularly complicated at small radii (< 10 AU). The emissivity has two peaks going down 
into the disk (increasing vertical column at fixed radius) that result from a combination 
of effects: density and temperature variations, abundance changes and line-trapping. The 
declining temperature, the conversion of atomic oxygen into molecules (CO, H 2 and O2), 
and trapping at large columns all tend to decrease the emissivity. These effects are illustrated 
by Figure dni which shows the specific emissivity (per O atom). For small column densities, 
it is constant when the population is thermalized (at radii R < 10 AU, where the critical 
densities are reached), or it increases toward a maximum when the critical densities are only 
attained at larger vertical column densities. Then it decreases again with vertical column 
density, except at the smallest radii (< 1 AU), where the rise in temperature due to viscous 
heating close to the mid-plane has the opposite effect. The occurrence of the two peaks in 
emissivity at intermediate columns is the result of the rapid increase of volumetric density 
with increasing depth. Significant emission occurs beyond columns of Ah = 10 22 cm -2 as a 
consequence of the fact that the excitation energy is relatively low (« 230 K) and contributes 
~ 10% of the total line flux. 

The l D 2 - 3 Pj 01 6300/6363 A and the % - 1 D 2 5577 A forbidden transitions are 
optically thin and sub-th ermally excited. In addition to the electronic rate coefficients 



( jZatsarinny fc Tayalll2003l ). t he collisional de-exc itation rate of the 1 D 2 level by atomic hy- 



drogen has been calculated by lKrems et al.l (120061 ). It varies little between 500-6000 K, where 



it has a mean value of 8 x 10~ 13 cm 3 s~ x to within 25%. At 4000 K, the collisional electron 
de-excitation rate is 3 x 10~ 9 cm 3 s _1 for the l D 2 level, so electronic collisions dominate over 
those with atomic hydrogen for x e > 2.5 x 10~ 4 . If we consider only electronic collisions, 
then an upper limit to the critical density for the 6300/6363 A transitions IS Tip ~ 10 cm , 
while that for the 5577 A transition is n e ~ 10 8 cm -3 . In the absence of calculations for the 
1 S'o level, we have adopted the value, 10~ 12 cm 3 s" 1 , for the H rate coefficient. Thus the 6300 
and 6363 A lines are populated near-thermally, and we can use the two-level formulae given 
in equations (3-2) and (3-3) of GNI07 that allow for the effects of finite values of the critical 
density. The emissivity for both lines is then given by equation [U The 1 S'o uppermost level 
is so far from being thermalized that we can calculate the flux of the 5577 A line by using 
the extreme low-density collisional-excitation limit, in which every excitation leads to the 
production of a photon. 

Figure [IT] shows that the O I forbidden lines trace the upper atmosphere of the disk, 
where the atomic oxygen abundance is constant. The OI 6300/6363 A lines are generated 
within 25 AU, a somewhat smaller range than the neon fine-structure lines, while the 5577 A 
line (not shown) originates from radii inside 15 AU. These calculations of forbidden line 
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fluxes necessarily ignore an effect of the UV photodissocation of OH, which leads to significant 
branc hing to the exited 1 D?, level o f 0, producing fluorescent emission of the 01 6300/6363 A 
lines ( IStorzer fc Hollenbachl |2000| ) . Because significant OH abundances only occur at large 
depths in our model, this effect might become important at large radii due to OH photo- 
dissociation by the stellar or interstellar /intracluster FUV radiation field. 



4.3. Sulfur 

We focus on the mid-IR fine-structure lines of SI and the 6718/6732 A forbidden tran- 
sitions of S II; the fine-structure lines of S III are likely to be too weak to be observable at 
the present time. In accord with Figure [HJ the calculation of the S I fine-structure emis- 
sion is similar to that just described for O I. The big difference is that little if anything is 
known about the excitation of these lines by collision partners other than electrons. Thus 
the present flux estimates are lower limits. The calculations take line trapping into account 
and solve the three-level fine-structure population problem exactly, including the the weak 
J = — 2 quadrupolar transition. 

Figure [T2l shows the spatial distribution of the S I 25.25 /im fine-structure line emissivity 
for the reference model. The appearance of the 56.33 fim emissivity (not shown) is similar, 
but it is concentrated toward smaller radii and higher altitudes. As in the case of O I, the 
variation of the fine-structure emission with depth at fixed radius is non-monotonic, and 
again there are several factors that contribute to this behavior. First, there are the chemical 
changes already mentioned, where sulfur changes from S + to S at high altitudes and from S 
to sulfur molecules at low altitudes. Thus the emission peaks at intermediate altitudes, at 
a vertical column of A H ~ 10 21 cm -2 at R = 0.25 AU and near A H ~ 10 19 cm -2 at R = 
25 AU. Second, the increase of volumetric density with depth also promotes fine-structure 
emission, whereas the overall temperature decrease has the opposite effect. Figure [TBI shows 
the variation of the specific emissivity with vertical column at fixed radius. For radii in the 
range 1-25 AU, the emissivity per S atom of the J = 1 — level changes from sub-thermal to 
thermal as the density increases with increasing depth. At somewhat larger column densities, 
the specific emissivity decreases due to decreases in both temperature and electron density. 

The transition of atomic sulfur into molecules such as SO and SO2 has little effect on 
the absolute level of the S I fine-structure emission. This is due to the fact that electronic 
excitation has already become ineffective for Ah < 10 22 cm~ 2 before molecule formation 
occurs. In order to gauge the importance of H atom collisions, we carried out a calculation 
where we assumed that the SI fine-structure de-excitation rates are about the same as for 
01, specifically jfe(H) = 5.0 x 10~ n T 0A for all transitions. The SI 25.25 and 56.33 /xm 
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emission is then increased by factors of 10 or more over that for electron collisions only. This 
serves to indicate the importance of extending the calculation of H excitation rates to heavy 
atoms and ions such as Ne II and S I. 

We estimate the emission of the S II forbidden lines in the optically thin approximation. 
We concentrate on the 2 -05/2,3/2 - 4 Sz/ 9 . transitio ns with wavelengths near 6700 A. Using 



the electron collisions strengths from iTayall (119971 ). we find electron critical densities, n cr ~ 
10 4 — 10 s cm -3 s -1 , and conclude that the lines are almost thermalized. We also considered 
the possible role of the H + + S — > S + + H charge exchange, which mainly branches to 
the upper 2 Pj level of S + , and then decays with the emission of photons with wavelengths 
near 4070 A, 6700 A, and 1 /mi. We find that this process does not increase the 6718 and 
6732 A line emissivities, probably because the rate coefficient for charge exchange is small 
(~ 5 x 10~ 12 cm 3 s _1 ) and because the high density tends to a thermalize the level population. 
In Figure [T4"t we show the emissivity of S II 6718 A line (the 6733 A line is very similar). These 
lines trace the warm upper layers of the disk. Significant emission arises only from regions 
with temperatures in the range 2000-5000 K. 



4.4. Carbon 

The energy levels of the ground configuration of C I are similar to those of O I and 
S I shown in Figure EJ except that the scale of the energy level separations is reduced 
and the ordering of the total angular momentum quantum number J in the ground state 
triplet is inverted. The upper fine-structure levels have energies E 2 /k = 62.5 K and E\/h = 
23.62 K and give rise to far-infrared lines with wavelengths A(l — 0) = 609. 135 /im and 
A(2 — 1) = 370.144 /mi. In addition to these lines, we calculate the emissivity of the forbidden 
lines at 9827 and 9853 A that emanate from the first 1 D 2 level of C I above ground. For 
CII, we only consider the fine-structure doublet, with excitation energy and wavelength, 
A£(3/2 - l/2)/k = 91. 2K and A(3/2 - 1/2) = 157.7 /mi, and ignore the 2300A lines that 
arise from the next 4 P level. 

The A values of the CI and CII fine-structure transitions are very small (~ 10 _7 s _1 
for CI and 2.3 x 10~ 6 s _1 for CII), which implies that the critical densities for all of these 
transitions are low (< 10 4 cm -3 , based on the collisional rate coefficients referenced in Table 
2). Thus we can estimate the emissivities assuming that the fine-structure levels are in 
thermal equilibrium (but we do include line-trapping). The integrated fluxes for a nominal 
distance of 140 pc are given in Table 3. The flux of the CII 158 /im line is very small, but 
the 369 /fm line of CI appears to be in the observable range, as may even be the case for the 
609 /im line. 
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We estimate the the CI 9827 and 9853 A forbidden line emission in the optically thin 
approximation. Using Table 2, the critical density of the 1 D 2 level at 4000 K is 1.96 x 
10 4 cm -3 s" 1 , and the level is close to being thermalized. The density-dependent correction 
for sub-thermal excitation is made in the same way as for O I, as described in Section 4.2. The 
spatial distribution of the C I 9827 A emission is shown in Fig. [151 The emission is peaked at 
small vertical column densities, much like the O I 6300 A emission in Fig. [Til The integrated 
line flux given in Table 3 is somewhat smaller than the flux of the O I 6300 A line, but still in 
the potentially observable range. Of course all of these statements about the detectability of 
neutral carbon lines stand to be corrected (and reduced) when interstellar/intracluster UV 
radiation is included in the ionization theory. 



5. Discussion 

In this section we discuss the most important results of this paper in the context of both 
existing and prospective observations. We treat high and low excitation lines separately, since 
they trace different parts of the disk. We first summarize in Table 3 the line fluxes for five 
values of the stellar X-ray luminosity that range from Lx = 2 x 10 29 — 2 x 10 31 erg s _1 . To 
ensure accurate results for values of Lx larger than standard, we have calculated the disk 
properties and emissivities out to 100 AU. The calculations pertain to the case where X-ray 
heating of the disk atmosphere dominates mechanical heating (a^ = 0.01 in the notation 
of GNI04). A s discussed in S e ction 2, the underlying density model is the generic T Tauri 



disk model of iD'Alessio et al.l (119991 ); the parameters for the model are given in the caption 
of Fig. 1. The "reference" model that we have used throughout has an X-ray luminosity of 
Lx = 2 x 10 30 erg s _1 ; the fluxes for this case are given in the middle column of the Table. 

The values shown in the other columns illustrate how sensitive the modeling results are 
to the choice of X-ray luminosity. Table 3 encompasses a range of 100 in X-ray luminosity, 
which is roughly the spread seen by the COUP and XEST projects for T Tauri stars in the 
Orion Nebula and Taurus-Aurigae clusters of the same mass and age. Were we to choose 



the median Lx for Taurus-Aurigae (jTelleschi et al.ll2007l ) for the reference model, the second 



column of fluxes would be appropriate. According to the discussion in the previous sections, 
all of the flux estimates are themselves uncertain to varying degrees due to uncertainties 
in the rate coefficients and other limitations of the model. In particular, the fluxes for the 
neon and sulfur fine-structure lines are lower limits because of the omission of collisional 
excitation by H atoms. 



The sensitivity of the integrated line fluxes to the X-ray luminosity can vary from one 
line to another, since they may originate in different parts of the disk with different physical 
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properties. We illustrate this in Fig. [T6], which shows how some of the line fluxes vary 
with Lx-i normalized to the standard model. The CI 369 /im line luminosity depends very 
weakly on Lx, since most of the emission is produced at intermediate column densities 
(iVn ~ 10 22 cm -2 ) where the X-ray flux is strongly attenuated. The 01 63 /im luminosity 
shows a somewhat steeper dependence, because it is produced both at high altitudes, where 
the X-ray flux is large, and close to the midplane, where the X-ray flux is small (and the gas 
and dust temperatures are almost the same). The optically thin O I 5577 A emission displays 
the strongest dependence on L x because, with its relatively high excitation temperature, it 
is only produced in the warmest part of the disk fully exposed to stellar X-rays. The other 
forbidden lines and the SI fine-structure lines (not shown), behave very much like the Nell 
fine-structure line. The almost linear relation between Nell line flux and Lx arises from 
the dependence of the emissivity on n\ (c.f. Eq. 4-1 of GNI07). In its present form, Fig. [TBI 
is not meant to suggest a general empirical correlation of line fluxes with X-ray luminosity, 
simply because the stellar and disk properties of the reference model have been held fixed 
in the calculations for this figure. 

5.1. High excitation lines 

5.1.1. Neon Fine- Structure Lines 

The neon fine-structure line fluxes (Table [3]) are close to but more accurate than GNI07. 
The Ne III flux is 40% larger since the calculations extend to smaller radii (R = 0.25 AU). The 
Nell 12.81 /xm line emission arises from the top of the atmosphere, where the temperature 
and electron density are high. We find significant contributions to the integrated emission 
out to R = 25 AU. The emissivity per unit radius is determined by the column density in 
the upper level (c.f. Eq. 3-7 of GNI07). This quantity is plotted in the upper panel of Fig. [T7] 
against the radius R. In contrast to Fig. 4 of GNI07, there is no dip near 5AU, thanks 
to more accurate and closer spaced calculations. When we assume Keplerian rotation to 
convert from column density to a velocity distribution function P(v) (c.f. Eq. 3-8 of GNI07), 
the peak in P(v) occurs at v — 0.25t> (1 AU), which corresponds to a peak in the emissivity 
per unit area at 16 AU. The bottom panel of Fig. [T71 shows that P(v) has a long tail extending 
to high velocity. The ratio of the Nelll 15.55 /xm to the Nell 12.81 /xm flux is of order 0.1. 
GNI07 considered the detection of the Nelll line to be a near definitive diagnostic of X- 
ray ionization of neon. However, using a disk model dominated by stellar EUV radiation, 
Hollenbach & Gorti (2007, private communication), find a range of possible Nelll/Nell ratios 
(0.1-6) depending on the the nature of the EUV spectrum. Thus it is not immediately clear 
whether the Nelll/Nell fine-structure line ratio can serve as a discriminant between EUV 
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and X-ray irradiation. 



The Nell 12.81 /im line has now been detected in T Tauri dis ks in a bout 20 cases from 
the ground an d from space (by the IRS on Spitzer. Pascucci et al. 120071; Lah uis et al. 120071 ; 
Espaillat et al. 120071 : and MICHELLE on Gemini North: Herczeg et al. 120071 ) . The detected 



flux levels are of order 10 -14 erg cm -2 s _1 ; they are within a factor of few agreement with 
our model calculations. The weaker 15.55 /im of NeHI line has been tentatively detected 
in only one case (Sz 102, Lahuis et al. 2007). Based on the 3a detection of the Nelll line, 
the Nelll/Nell ratio is ~ 0.06, which might be compared with a ratio ~ 0.1 derived from 
Table 3. Perhaps what is most interesting about this T Tauri star is that it has the second 
brightest Nell line detected in the Lahuis et al. survey. Its flux is an order of magnitude 
larger than the value for our reference model, after correction for the distance of 200 pc. 
Detailed modeling of th e disk around this st ar would be of considerable interest. For the 
case of nearby TW Hya, iHerczeg et al.1 (120071 ) measure the luminosity for the Nell 12.81 /im 
line to be 4.8 x 10~ 6 L Q , which is 20% larger than the result for our reference model in Table 
3. 

From a small number of detections of the Nell 12.81 fim li ne the Spitzer Lagacy pro- 
gram, "Formation and Evolution of Planetary Systems" (FEPS). lPascucci et al.1 (120071 ) have 
suggested that a linear correlation may exist between the Ne II and X-ray luminosities. To- 
gether with their detection of the Ne II line in CS Cha and that for TW Hya, based on 
data obtained by lUchida et al.l (120041 ). lEspaillat et al.l (120071) argue against suc h a corre- 
lati on. The range of Lv in these papers is small, 1.5 for iPascucci et al.l (120071 ) and 2.25 
for lEspaillat et al.l (120071 ) , the same order of magnitude as the observed variability in the 
X-ray luminosity. Aside from the intrinsic difficulty of deducing a meaningful correlation on 
the basis of a small dynamical range in YSO X-ray luminosity, it is important to keep in 
mind that the line fluxes depend on many other stellar and disk properties, e.g., the stellar 
temperature, mass, and age and the disk mass, accretion rate, and the amount of dust grain 
growth and settling. The theoretical calculations shown in Fig. [16] may not be relevant in 
this context because they assume that the disk structure is fixed while Lx is varied. To 
pursue the question of the existence of a correlation with Lx, T Tauri stars with similar 
properties, both stellar and disk, should be considered. A more practical approach might be 
to focus in depth on a few specific cases with well observed properties and to use observed 
X-ray spectra. 



Using the MICHELLE spectrometer at Gemini North, IHerczeg et al.l (120071 ) have ob- 
tained a spectrally resolved Ne II line profile for TW Hya, which has a nearly face-on disk. 
The line is centered at the stellar radial velocity, consistent with a disk origin, but the line 
is broad with FWHM of ~ 22 km s~ x . The rotational broadening in our models (Fig. 5 of 



-18- 



GNI07 and Fig. [T7] of this paper), when corrected for inclination, cannot explain the ob- 
served line width, even though the present calculations extend to radii as small as 0.25 AU. 
Although our calculations agree with the measured Nell luminosity, the rotational broad- 
ening in the model is small because it arises from relatively large radii, as shown in Fig. [T71 
As discussed by Herczeg et al. (2007), two interesting explanations of the line width seen 
TW Hya that deserve consideration are transonic turbu lence and photo-evaporative outflows 
from the disk (jHollenbach et al.ll2000l ; iFont et al.l 120041 ) . These possibilities raise interesting 
challenges for future research. If the hole is produced by photo-evaporation, then the wind 
itself will generate line emission. Furthermore, the i nner rim of the hole will be irr adiated 
more or less directly by the star, as in the model of IChiang fc Murray- Clayl (120071 ). again 
producing characteristic line emission. 



5.1.2. Sulfur Fine- Structure Lines 

We consider the S I fine-structure lines to be "high-excitation" because the associated 
excitation energies, 476 K and 825 K, are larger than the temperatures characteristic of the 
outer disk, as are the fine-structure levels of the neon ions. As discussed in Section 3, 
the X-ray ionization theory for neon and sulfur, where charge exchange plays an important 
role, are also similar. Reference to Table 3 shows that the integrated S I fine-structure 
emission for our reference model is ~ 20 smaller than for the Ne II line. Fig. [T2] shows that 
the strongest emission of the SI 25.25 /im line comes from a relatively thin layer between 
Nr = 10 19 — 10 21 cm~ 2 for R < 10 AU. This situation and the attendant weakness of the 
25.25 /im line occur because sulfur is in the form of S + at higher altitudes and electronic 
excitation becomes weak at lower altitudes (c.f. Figs. [5] and E]). 



The S I 25.25 /im line was not detected by the Spitzer c2d and FEPS teams. iPascucci et al. 



(120071 ) put upper limits on the residual gas in the disks around their sur yey of older T Tauri 



stars w ith ages between 5 and 100 Myr using the model calculations of iGorti fc Hollenbach 



(120041 ) . From Table 3, which is based on pure electronic excitation, the S I fine-structure lines 
from younger T Tauri disks would be similarly difficult to detect. However, if the atomic H 
collision rates for S I are similar to those for O I, the predicted fluxes would approach the 
observable range. It is important to recall several other caveats that pertain to our calcula- 
tions: the uncertainties in sulfur chemistry, the total abundance of gaseous sulfur in disks, 
and the poorly known charge-exchange rate coefficients. 
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5.1.3. Forbidden Lines 



Table [3] gives the fluxes for our reference model of a sample of forbidden lines, as cal- 
culated in Section 4: 01 6300/6363 A, 01 5577 A, S II 6718/6731 A; CI 9827/9853 A. We 
have identified them as potential diagnostics of the warm and highly-ionized gas produced 
by X-ray irradiation in the upper atmosphere of the inner disk. These lines are generated 
within the same radial distance range as the Nell 12.81 and 15.55 /im lines, but at somewhat 
higher elevations. In principle, they can provide complementary checks on our picture of X- 
ray irradiation. As discussed in Section 4, the accuracy of the calculations varies somewhat 
from line to line, depending on the reliability of the atomic data and on the completeness of 
the model itself. For example, we have ignored the role of stellar or interstellar/intracluster 
FUV radiation, which could enhance the ionization of carbon at high altitudes and thereby 
reduce the amount of neutral carbon and the strength of the C I forbidden lines. Similarly, 
the SII 6718 A line strength is affected by the sulfur atoms being converted into molecules. 
With this proviso, we can see from Table [3] that many of the forbidden lines generated in the 
upper disk atmosphere are observable in nearby star-forming regions such as Taurus- Aurigae 
in the sense that the predicted fluxes are > 10~ 15 erg cm" 2 s" 1 , notably 01 6300/6363 A, 
CI 9827/9853 A, and perhaps SI 6718/6733 A. 



Many of these lines have been measured in T Tauri stars. The most dramatic manifes- 
tation of t he for bidden lines is their tracing of high- velocity outflows or jets from YSOs (e.g., 
Ray et al. 120071) . But these flows have a low- velo ci ty as well as a high-velocity component 
flKwan fc Tademarulfl98~8l . Il995l : lliirth et alill997T l. iKwanl Jl997T l suggested that the low ve- 



locity component (LVC) originates in the inner regions of T Tauri disks (R < 2 AIT 
possible sources for the origin of the LVC is a MH D disk wind (e.g., Pudritz et al 



j Ot her 



a photoevaporative outflow (e.g., Hollenbach et al. l2000l ; Font et al. 120041 ) 



20071 ) or 



Hartigan et all (119951 ) have obtained extensive data on forbidden line emission from 
T Tauri stars. The LVC is ubiquitous, with a typical line luminosity of ~ 1O~ 4 L , an 
order of magnitude larger than the luminosity of our reference model, 8.0 x 10~ 6 L Q . They 
also find OI 5577 A to OI 6300 A line ratios in the range ~ 0.2 — 0.5; our model predicts 
~ 0.1. These numbers may signify the difficulty of applying our X-ray irradiated disk model 
to observations such as the optical forbidden lines that are sensitive to the other dynamic 
entities that are operative in accreting young stars, e.g., the various flows that arise near 
the inner edge of the disk. Hartigan et al. define the LVC to include velocities as high as 60 
km s _1 , thereby increasing the possibility of including emission from such flows. A cleaner 



3 Kwan's theory of warm disk coronae has similarities with our X-ray irradiated disk atmospheres, but 
there are significant differences in the underlying ionization physics. 
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comparison case for our model might be provided by an X-ray bright star with a relatively 
low accretion rate, where the emission from outflows would be reduced. 

Because TW Hya has a low mass-loss r ate an d a small accretion rate, it would be a good 



case for detailed modeling. iHerczeg et al.l (120071 ) detected the 01 6300/6363 Alines in TW 
Hya at moderately high spectral resolution. The luminosity of the 6300 Aline, 7.0 x 10~ 6 L Q , 
is close to the result given in Table 3, 8.0 x 1O _6 L , recalling the good agreement for the 
luminosity of the Nell 12.81 //m line. The 6300 A line is centered on the stellar velocity, 
and it has a significantly narrower width than the Nell line, but one that is still too large 
to be explained by pure rotational broadening. Again, interesting possible explanations to 
consider in future modeling are a turbulent atmosphere and photo-evaporation. 



5.2. Low excitation lines 



We next consider the low-excitation lines that arise from the fine-structure levels of O I 
(excitation energies 228 K and 327 K), CI (excitation energies 23.6 K and 62.5 K) and CII 
(excitation energy 91.2 K); they give rise to far-infrared emission at 63 and 145 /im (OI), 
609 and 370 /im (CI), and 158 /xm (CII). Table 3 shows that the OI fine-structure lines are 
potential diagnostics of both the moderately warm, ionized gas produced by X-rays, and 
the cooler gas at large perpendicular column densities. As shown in Fig. [9J O I 63 /im line 
emission extends down toward the mid-plane region, especially at small radial distances, 
due to the incomplete conversion of atomic oxygen into molecules and the relatively low 
excitation temperature of the J = 1 — 2 transition. 

The 63 /im line has bee n obse rved around T Tauri stars with the K AO (e.g., Cohen 



et al. Il988l ; Ceccarelli et al. I1997T) and ISO (e.g, Spinogolio et al. |2000| ; Creech-Eakman 
et al. l2002i ; Liseau et al. 120061 ). These observations were made with large beams (tens of 
arc seconds) and with moderate spatial resolution (50-300 km s _1 ). A typical flux for a T 
Tauri star in Taurus- Aurigae is quite large, ~ 10~ 12 — 10 _11 erg cm -2 s -1 , and may include 
emission from other circumstellar material and not just disks. On the other hand, the flux 
calculations in Table 3 probably under estimate the flux, eve n though they integrate out to 
100 AU, especially for the carbon lines. IJonkheid et al.l (120041 ) have applied a code developed 
for (UV) photon-dominated regions to the surface layers of a flaring disk, and find that the 
fluxes of the O I, CI, and C II fine-structure lines are all of comparable brightness, and not 
that much weaker than the CO(J =1 — 0) line. Photo-ionization of neutral carbon by 
external FUV radiation can enhance the abundance of C + and increase the flux of the C II 
158 /im line. Our model would have to be extended to include UV irradiation in order to get 
a fuller understanding of these low-excitation diagnostic lines. 
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Our estimate of the O I 63 jjm line flux provides a good estimate of the emission of 
this line from the inner disk. Since its detection requires airborne or space observations, 
the identification with inner disk emission requires resolved line profiles and the assumption 
of Keplerian disk rotation. Such measurements could be made with heterodyne detectors 
aboard SOFIA or with Herschel (PACS). 



6. Summary 

We have used a simplified thermal-chemical and structural model to explore the diag- 
nostics of the atomic gas in protoplanetary disks irradiated by stellar X-rays. Our focus has 
been on moderate disk radii out to roughly 25 AU, where stellar X-rays play an important 
role in ionizing and heating the upper layers of the disk. We find several diagnostics for this 
region, all of which arise from relatively high-excitation levels ~ 1000 K: Nell and Nelll fine- 
structure lines and O I, S II, and C I forbidden lines. Some low-excitation fine-structure lines 
are also of interest in the context of X-ray irradiation, especially the O I fine-structure lines, 
which generates strong disk emission over a substantial range of vertical column density, 
extending close to the mid-plane of the disk. 

The atomic lines probe different depths of the disk at moderate radii. The forbidden 
transitions are generated at the top of the atmosphere, where the temperature and ionization 
levels are the highest. The mid-infrared fine-structure lines probe greater depths, where the 
temperature and ionization level have decreased. Neon is the cleanest case, since it does not 
easily form molecules, and its fine-structure lines originate from intermediate as well as top 
layers of the disk. On the other hand, the SI fine-structure lines are produced mainly in 
intermediate layers, since atomic sulfur has a low abundance in the top layers due to X-ray 
ionization and also near the mid plane, due to molecule formation. By contrast, the 01 
fine-structure lines are formed at almost all depths, including close to the mid plane, where 
the density is high and the upper levels of the transitions are more commensurate with the 
temperatures there. 

The X-ray generated atomic lines discussed here are produced over a range of radii, 
extending from the inner radius of the disk out to 20-30 AU. To resolve these regions in 
nearby star-forming clusters at a nominal distance of 140 pc, requires a spatial resolution 
of ~ 0.1 — 0.2", as well as resolved line profiles. For example, the Nell 12.81 /im line 
should be detectable with TEXES or comparable spectrometers on large ground-based tele- 
scopes. A promising step in this direction is its detection with MICHELLE on Gemini 



North (IHerczeg et al.l 120071 ). The far-infrared lines of CI, CII, and OI will be detectable 



with instruments on board SOFIA and Herschel. 
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In presenting these results for a simplified T Tauri disk model irradiated by stellar X- 
rays, we have tried to mention and discuss its limitations. To summarize, there are, first of 
all, many uncertainties in the underlying atomic physics, such as missing rate coefficients for 
charge exchange with atomic hydrogen and for collisional excitation by abundant particles 
other than electrons. These uncertainties of course hold for all thermal-chemical/excitation 
models. Another limitation is the reliance on one type of disk, represented by the smooth 
density distribution of the generic T Tauri model of D'Alessio et al. (1999). Removing this 
restriction is one of our immediate goals. It will enable us to discuss disks of different ages 
and disks with non-monotonic density distributions such as holes, gaps, and rims. The 
improvements that are needed to treat situations such as disks with holes require a higher 
level model than used here. More generally, a 3-dimensional treatment is needed, and this 
dictates a greatly improved treatment of the radiation transfer, especially for the cooling 
lines and the external UV irradiation. 

Focusing on X-rays to the exclusion of UV irradiation is another limitation that needs to 
be dealt with in future modeling. It is very likely that the well-documented X-rays produced 
by essentially all YSOs are the dominant external radiation source for radii less than 25 AU. 
At larger radii, however, ambient UV, from the general interstellar medium or from the host 
star cluster members, can affect disk properties. At smaller radii, UV radiation originating 
close to the star may also play a role, although it is more easily attenuated than moderately 
hard keV X-rays. We have ignored stellar UV radiation, which is poorly determined, in favor 
of the well-measured X-rays. 

Granted that the X-rays likely play the dominant role within 25 AU or so, the relatively 
high-excitation transitions that they generate may be produced by some of the other dynam- 
ical components of the star-formation system that lie near and even inside the inner edge of 
the disk. The most obvious examples are the forbidden optical transitions produced by jets 
and accretion streams that are the signatures of actively accreting YSOs. In this situation, 
the high-excitation line strengths that we calculate for disks may only represent lower limits 
to those measured in spatially unresolved observations. A similar conclusion pertains to the 
low-excitation lines, where significant emission can arise from large disk radii where external 
UV can dominate. Of course the contributions of the various parts of the star formation 
complex can in principle be disentangled by observations with the appropriate spatial and 
spectral resolution, together with more complex models. The present results demonstrate 
the importance of X-ray irradiation in the ongoing process of elucidating the nature of the 
gas in protoplanetary disks. 
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Fig. 1. — Gas temperature structure of the reference model described in the t ext w ith color 
code on the right. The density structure is based on the iD'Alessio et al.l (119991 ) model, 
defined by the parameters: M* = O.5M , R = 2R & , T* = 4000 K, M = 10" 8 il^yr" 1 . The 
reference X-ray luminosity is Lx = 2 x 10 30 erg s _1 , and the temperature of the thermal 
X-ray spectrum is kTx = IkeV. The white curves are contours of fixed vertical column 
density Ah in units of cm -2 . The curve for JVh = 10 21 cm -2 corresponds roughly to two 
scale heights as defined near the mid-plane. The integrated surface density varies roughly 
as 1/R for R > 1 AU and somewhat less rapidly for R < 1 AU. 



Electron abundance 




Fig. 2. — Spatial distribution of the electron fraction plotted in the same way as Fig. 1, 
except for the color-coded units. 



-29- 



Ne + fractional abundance 




10 20 30 40 



R (AU) 




R (AU) 

Fig. 3. — Spatial distribution of the Ne + and Ne 2+ abundances relative to the total abun- 
dance of neon, plotted in the same way as Figure 1. 




Fig. 4. — Schematic diagram of sulfur chemistry. Charge exchange of sulfur ions with 
hydrogen atoms is critical for the ionization structure at high altitudes, whereas molecule 
formation controls the abundance of atomic sulfur near the mid-plane. 
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Fig. 5. — Fractional sulfur abundances vs. vertical column density measured from the top 
of disk at a radial distance of 20 AU from the star. SO2 becomes the dominant species at 
large columns. 
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Fig. 6. — Spatial distribution of normalized S and S + abundances, plotted in the manner of 
Figure 1. 
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Fig. 7. — The spatial variation of the Nell 12.81 /xm fine structure emissivity, defined by 
equation (1), in units of erg cm _3 s _1 . 



-32 - 



Ol 



SI 



48620' 



22830 



{ 



326- 
227- 



2972/ 
2958A 



- 1 s. 



5577A 



31912 



_1 D 



6363/6300A 



tl45.5un 



L63.18jjm 



13292 



#■825 
J 570 

I 



4591 A 


7727A 






1.13/1.08|jm 






t56.33|jm 




| 25.25|jm 



Fig. 8. — Energy level diagrams for O and S showing the low-lying levels that generate the 
fine-structure and forbidden lines. Note that the energy scale is in Kelvin. 
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Fig. 9. — Spatial distribution of the O I 63 /zm fine-structure line emissivity, defined by 
equation (1), in units of erg cm~ 3 s _1 . 
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Fig. 10. — The specific 01 63 /im emissivity in units of erg s 1 atom , vs. perpendicular 
column density at various radii. 



I 6300 A emissivity 




R (AU) 

Fig. 11. — Spatial distribution of the 01 6300 A forbidden line emissivity, in units of erg 
cm s . 
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Fig. 12. — Spatial distribution of the SI 25.55 /zm fine-structure line emission in units of erg 
cm -3 s _1 . The emission peaks at intermediate depths, as discussed in the text. 
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Fig. 13. — The spatial distribution of the specific emissivity in units of erg cm 3 s 1 atom 
of the SI 25.55 /xm fine-structure line vs. perpendicular column density at several radii. 
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Fig. 14. — Spatial distribution of the S II 6718 A forbidden line emission in units of erg cm 
s _1 . Significant emission only occurs in the warmest part of the disk atmosphere. 
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Fig. 15. — Spatial distribution of the CI forbidden line emission in units of erg cm 3 s 
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Fig. 16. — The variation of selected line fluxes with X-ray luminosity, all normalized to the 
values for the reference model. 



-37- 



100.000 




w=v/v(l AU) 

Fig. 17. — Column densities of excited Ne ions in units of 10 14 cm~ 2 plotted vs. radial dis- 
tance in AU (left), and the rotational velocity distribution function P(v) plotted vs. velocity 
normalized to its value at 1 AU (right). 



Table 1: Selected Reaction Rate coef 


icients for Sulfur 


React ants 


Products 


Symbol 


Rate coefficient 


Reference 


S+ + e 
S 2+ + e 
S 3 + + e 


S + hv 
S + + hv 
S 2+ + hv 


a 2 
a 3 


4.65 x iCr 13 T 4 -°- 63 
3.37 x i(T 12 T 4 - a50 
5.23 x i(r 12 T 4 - a50 


Aldrovandi & Peauimot fl973): Gould ( 1978) 
Nahar & Pradhan ( 1995) 
Nahar & Pradhan ( 1995) 


H+ + S 
S+ + H 
S 2 + + H 
S 3 + + H 


S+ + H 
H+ + S 
H+ + S+ 
H+ + S 2+ 


ki 
h 

h 


5 x iO" 12 * 
negligible 

1Q -14 

3 x iO" 9 


Zhao et al. (2005) 
Zhao et al. (2005) 

Butler & Dakarno (1980): Christensen & Watson (1981) 
Butler & Dakarno (1980): Christensen & Watson (1981) 


S 2 + + H 2 
S 2 + + H 2 


S etc. 
S+ etc. 


BK 

(1 - B)K 


2 x fO" 9 B 
2 x iO- 9 (f-B) 


Chen et al. (2003) 
Chen et al. (2003) 


S+ X-ray 
S+ X-ray 
S++ X-ray 


S+ + e 
S 2+ + 2e 
S 2 + + e 


Csec(S) 

C dir (S) 
C(S) 


5xC 
20 xC 
25 x C 


This work. 
This work. 
This work. 



* Approximate value valid from 500-5,000 K 



Table 2: References for excitation coefficients* 



Ion 


references 


Ne II fine structure 
Ne III fine structure 
1 fine structure 

1 forbidden 
S I fine structure 
S II forbidden 
C I fine structure 

CI forbidden 

C II fine structure 


Griffin et al. (2001) (e) 
Butler & ZeiDDen (1994) (e) 

Launav & Roueff (1977a): Abrahamsson et al. (2007) (R): Jaauet et al. (1992) (Ho): 
Monteiro & Flower (19871 (Re): Chambaud et al. (19801 (d): Peauimot (1990) (e) 
Krems et al. (2006) (H): Zatsarinnv & Taval (2003) (e) 
Tayal (2004) (e) 
Tayal (1997) (e) 

Launav & Roueff (1977a); Abrahamsson et al. (2007) (H); 
Monteiro & Flower (1987) (He. Ho) 

Peauienot & Aldrovandi (1976); Zatsarinnv et al. (2005) (e) 

Barinovs et al. (2005) (H): Launav k Roueff (1977b) (Ho): Wilson & Bell (2002) (e) 



*The collision partners are specified inside the brackets. 



Table 3: Unresolved line fluxes in ergs cm 



L x (erg s" 1 ) 


2.0 x 10 29 


6.0 x 10 29 


2.0 x 10 30 


6.0 x 10 30 


2.0 x 10 31 


Nell 12.82 jim 
Nelll 15.55 /an 


4.1 x 10~ 16 
3.0 x 10~ 17 


1.5 x 10~ 15 
1.4 x 10~ 16 


6.4 x 10" 15 
7.4 x 10" 16 


2.3 x 10" 14 
3.3 x lO" 15 


9.6 x 10- 14 

1.7 x 10~ 14 


T 63 nvn 

vy j. i (.111 

01 146 /xm 
01 6300 A** 
01 5577 A 


1.9 x 10~ 14 

9.6 x 10~ 16 
1.1 x 10~ 15 

2.7 x 10~ 17 


3.1 x 10~ 14 

1.2 x 10" 15 
3.6 x 10" 15 
1.4 x 10~ 16 


6.6 x 10 14 
1.8 x 10 15 
1.3 x 10 14 
7.3 x 10" 16 


1.6 x 10" 13 
3.8 x 10^ 15 
4.5 x 10~ 14 
2.8 x 10" 15 


4.3 x 10" 13 
1.1 x 10~ 14 
1.7 x 10~ 13 
1.1 x 10~ 14 


SI 25 25 urn 
SI 56.23 /xm 
SII 6718 A 
SII 6733 A 


7.0 x 10" 17 
6.4 x 10~ 18 

4.1 x 10~ 17 

9.2 x 10~ 18 


1.3 x 10~ 16 
9.0 x 10~ 18 
1.7 x 10~ 16 
3.9 x 10~ 17 


3.4 x 10 16 
1.7 x 10 17 
8.3 x 10 16 
2.0 x 10" 16 


9.5 x 10~ 16 
4.0 x 10~ 17 

3.5 x 10" 15 

8.6 x 10~ 16 


3.1 x 10" 15 

1.2 x 10~ 16 
1.6 x 10~ 14 
4.2 x 10~ 15 


CII 158 /xm 
CI 369 /xm 
CI 609 /xm 
CI 9827 A** 


1.2 x lO" 18 
9.8 x 10~ 16 
3.2 x 10~ 16 
3.5 x 10~ 16 


3.2 x 10~ 18 
1.2 x 10~ 15 
3.8 x 10~ 16 
1.4 x 10" 15 


9.8 x 10" 18 
1.6 x 10 15 
4.5 x 10 16 
6.3 x 10" 15 


2.8 x 10~ 17 
2.0 x 10~ 15 

5.4 x 10~ 16 

2.5 x 10~ 14 


9.0 x 10~ 17 
2.6 x 10~ 15 
6.4 x 10" 16 

1.1 x 10~ 13 



* A nominal distance of 140 pc has been assumed. Luminosities can be obtained by 
multiplying by 47r(140pc) 2 = 2.35 x 10 42 cm . The middle column refers to the 
reference model. 



** The 01 6363 A and CI 9853 A fluxes can be obtained by multiplying the fluxes of 
01 6300 A by 0.244 and CI 9827 A by 0.248, respectively. 



